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We use experiments, numerical simulations, and theoretical analysis to investigate the propaga- 
tion of highly nonlinear solitary waves in periodic arrangements of dimer (two-mass) and trimer 
I (three-mass) cell structures in one-dimensional granular lattices. To vary the composition of the 

, fundamental periodic units in the granular chains, we utilize beads of different materials (stainless 

' steel, bronze, glass, nylon, polytetrafluoroethylene, and rubber). This selection allows us to tailor 

^Sj ' the response of the system based on the masses, Poisson ratios, and elastic moduli of the compo- 

, , nents. For example, we examine dimer configurations with two types of heavy particles, two types 

^ ■ of light particles, and alternating light and heavy particles. Employing a model with Hertzian inter- 

actions between adjacent beads, we find very good agreement between experiments and numerical 
simulations. We find equally good agreement between these results and a theoretical analysis of the 
■ model in the long-wavelength regime that we derive for heterogeneous environments (dimer chains) 

^S) ' and general bead interactions. Our analysis encompasses previously-studied examples as special 

cases and also provides key insights on the influence of heterogeneous lattices on the properties 
(width and propagation speed) of the nonlinear wave solutions of this system. 
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I. INTRODUCTION 



Ever since the Fermi-Pasta-Ulam problem was first investigated half a century ago, nonlinear oscillator chains have 
J> ' received a great deal of attention [H, 0, H, 01 • Over the past few years, in particular, such chains of nonlinear oscillators 
OA , have proven to be important in numerous areas of physics-including Bose-Einstein condensation in optical lattices in 
atomic physics [1, d, 01 j coupled waveguide arrays and photorefractive crystals in nonlinear optics [1, Q , and DNA 
[ double-strand dynamics in biophysics [lo| . 

The role of "heterogeneous" versus "uniform" lattices and the interplay between nonlinearity and periodicity has 



' been among the key themes in investigations of lattice chains 11 1. Here we investigate this idea using one-dimensional 
(ID) lattices of granular materials, which consist of chains of interacting spherical particles that deform elastically 



when they collide. The highly nonlinear dynamic response of such lattices has been the subject of considerable 
attention from the scientific community [H, [H, [11 [11 [H, [13, [S [ll [13, [SI, [H, [H [H [H, [M, [13 . In contrast 



to the traditionally-studied nonlinear oscillator chains, these systems possess a "double" nonlinearity arising from 
the nonlinear contact interaction between the particles and a zero tensile response. This results in an asymmetric 
potential, which has in turn led to the observation of several new phenomena. In particular, some of the most 
d interesting dynamic properties of such highly nonlinear systems appear when the materials are under precompression 
[l2, 13, 14], at the interface between two different highly nonlinear systems [I^ [13, [Hlj [13 1 interface between 

linear and nonlinear structures (23j . This completely new type of wave dynamics opens the door for exciting new 
fundamental physical phenomena and has the potential to ultimately yield numerous new devices and materials. 
Because granular lattices can be created from numerous material types and sizes, their properties are extremely 
tunable [T^, [Tsl . [T3|. Furthermore, the addition of a static precompressive force can substantially vary the response 
of such systems and allows the selection of the wave's regime of propagation between highly nonlinear and weakly 
nonlinear dynamics [11]. Such tunability is valuable not only for studies of the basic physics of granular lattices but 
also in potential engineering applications. Proposed uses include shock and energy absorbing layers [TtI . [^T 
sound focusing devices (tunable acoustic lenses and delay lines), sound absorption layers, sound scramblers [1 
and more. 

In the present work, we focus on one of the fundamental physical properties of such chains. Namely, because they 
are highly nonlinear, granular lattices admit a novel type of wave solution whose qualitative properties differ markedly 
from those in weakly nonlinear systems. Indeed, it was the experimental realization of such waves and the th eory 
developed for uniform lattice systems of this type that has motivated the current broad interest in such settings [l2] ■ 
The remarkable property of these waves is that they essentially possess a support that consists of of just a few lattice 
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sites [12|, providing perhaps the closest experimentally tractable application of the notion of "compactons" [29| and 
establishing the potential for the design and creation of systems with unprecedented properties. 

Our concern in this paper, which provides a detailed and expanded description of work we reported in a recent 
letter [305, is to extend the established theory for uniform granular lattices to nonuniform ones. One way to do this 
is to study the effects of defects, such as inhomogeneities, particles with different masses, and so on. This has led to 
the observation of interestin g p hysical responses such as fragmentation, anomalous reflections, and energy trapping 
[TtI. [isl. [19I [20I [2ll. [2^. [2^. [2^.f25| . In the present paper, on the other hand, we examine the prototypical heterogeneities 
of granular lattice "dimers" and "trimers," consisting of chains of multiparticle cells composed, respectively, of two and 
three different types of beads. Such heterogeneities arise in a diverse array of physical settings, including ferroelectric 
perovskites [3l|,[3l| and polymers [11], optical waveguides [13, [IHl, and cantilever arrays [36l|. In our granular setting, 
we use a variety of soft, hard, heavy, and light materials to investigate the effects of different structural properties in 
the fundamental components of such systems. We also vary the number of beads of a given type in each cell in order 
to examine the effects of different unit cell sizes (i.e., different periodicities). 

More specifically, we investigate solitary wave propagation using experiments, numerical simulations, and theoret- 
ical analysis. We report very good agreement between experiments and numerics. For the case of dimer chains, we 
also construct a long-wavelength approximation to the nonlinear lattice model to obtain a quasi-continuum nonlinear 
partial differential equation (PDE) that provides an averaged description of the system. We obtain analytical expres- 
sions for wave solutions of this equation and find very good qualitative agreement between the width and propagation 
speed of these solutions with those obtained from experiments and numerical simulations. 

The rest of this paper is organized as follows. First, we present our experimental and numerical setups. We then 
consider chains of dimers and discuss our experimental and numerical results. We subsequently use a long-wavelength 
approximation to derive a nonlinear PDE describing the dimer setup in the case of general power-law interactions 
between materials and construct analytical expressions for the propagation speed, width, and functional form of its 
solitary wave solutions. We then compare this continuum theory to our experiments and numerical simulations of the 
discrete system. Finally, we discuss our experiments and numerical simulations for chains of trimers and summarize 
our results. 



II. EXPERIMENTAL SETUP 

The experimental dimer and trimer chains were composed of vertically aligned beads in a delrin guide that contained 
slots for sensor connections or in a guide composed of four vertical garolite rods arranged in a square lattice [see 
Fig. [Ija)]. Each 'Wi : N2 dimer" consisted of a variable number Ni g {1, . . . , 7} of one type of bead alternating with 
N2 e {1, . . . , 7} of a second type of bead in a periodic sequence. The 'Wi : N2 : N3 trimers" we studied are defined 
analogously. The predominant class of configurations we considered included A^i high- modulus, large mass stainless 
steel beads (non-magnetic, 316 type) and N2 low-modulus, small mass polytetrafluoroethylcnc (PTFE) or Neoprene 
rubber elastomer beads. We also examined steehbronze, PTFEiglass, and PTFE:nylon dimers and 1:1:1 trimers of 
steel:bronze:PTFE, steel:glass:nylon, and steel:PTFE:rubber. The diameter of all spheres was 4.76 mm. 

For each experiment, we connected four calibrated piezosensors to a Tektronix oscilloscope (TDK2024) to detect 
force-time curves. Three of the piezo-sensors were embedded inside particles in the chain and a fourth one was 
positioned at the bottom (i.e., at the wall). To fabricate the sensors, beads selected from the various materials were 
cut in two halves and slots for sensors and wires were carved within them. Vertically poled lead zirconate titanate 
piezo elements (square plates with 0.5 mm thicknesses and 3 mm sides), supplied by Piezo Systems, Inc. (RC=10'^ 
microseconds), were soldered to custom micro- miniature wires and glued between the two bead halves [see Fig.[ljb)]. 
We calibrated the setup using conservation of momentum. Waves were generated in the chains by dropping a striker 
from various heights. In most of our experiments, the striker consisted of a stainless steel bead; for some configurations, 
we also used a PTFE bead, a glass bead, and a PTFE:glass dimer. The material properties for the various beads are 
shown in Table [J 
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Material 


mass 


E 


V 


Steel 


0.45 g 


193 GPa 


0.3 


PTFE 


0.123 g 


1.46 GPa 


0.46 


Rubber 


0.08 g 


30 MPa 


0.49 


Bronze 


0.48 g 


76 GPa 


0.414 


Glass 


0.137 g 


62 GPa 


0.2 


Nylon 


0.0612 g 


3.55 GPa 


0.4 



TABLE I: Material properties (mass, elastic modulus i?, and Poisson ratio v) for stainless steel [37l. Is^. PTFE p^. l39l. l40l] . 
rubber (McMaster-Carr), bronze [U, |42|, glass [l^l, and nylon 44]. The value of the dynamic elastic modulus of the rubber 
beads was extrapolated from the experimental data. 



[a) 




FIG. 1; (Color online) (a) Experimental setup for a dimer chain consisting of a periodic array of cells with A''i consecutive 
beads of one material (e.g., stainless steel) and consecutive beads of another material (e.g., PTFE). (b) Schematic diagram 
of the composition of the sensors placed in the chain. 



III. NUMERICAL SIMULATIONS 

We model a chain of n spherical beads as a ID lattice with Hertzian interactions between beads [T2j : 

where j € {1, • • • , n}, is the coordinate of the center of the jth particle, bj = max{j/j_i — i/j, 0} for j e {2, . . . , n}, 
5i = 0, 5„+i = max{y„, 0}, g is the gravitational acceleration, Ej is the Young's (elastic) modulus of the jth bead, is 
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its Poisson ratio, rrij is its mass, and Rj is its radius. The particle j = represents the steel striker, and the (n + l)st 
particle represents the wall (i.e., an infinite-radius particle that cannot be displaced). Our numerical simulations 
incorporate the nonuniform gravitational preload due to the vertical orientation of the chains in experiments but do 
not take dissipation into account. 

The initial velocity of the striker is determined by experiments and all other particles start at rest in their equilibrium 
positions, which are are determined by solving a statics problem {yj ~ for all j). Each bead experiences a force 
from gravity (one can also include a constant precompression /-,-, but that is zero for the problem we study). One 
starts by considering the bottom of the chain with particle n against the wall. It experiences a force F\ from the top 
that is the sum of all the forces acting through the center of particle n and a force from the wall (particle n + 1) 
acting through the nonlinear spring (the Hertzian interaction). The expressions for F\ and i^2j 

n 

i=i 

F2 = m„A„,„+i(u„ - u„+i)^/^ , (2) 

must be equal at equilibrium. The displacement of the wall = is known, so Eq. can be solved for the 

equilibrium displacement of particle n in terms of known quantities: 




(3) 

One then repeats this analysis one particle at a time (from the bottom of the chain to the top) to obtain the other 
equilibrium positions. 



IV. CHAINS OF DIMERS 



Let us first discuss our results for chains of dimers. While we focus our presentation on steelrPTFE dimers, we also 
investigated steehrubber, steehbronze, PTFE:glass, and PTFEmylon dimers. Additionally, our presentation includes 
discussions of 1:1, A^i : 1, and 1 : Ni dimers. To compare numerical simulations with experiments, for which the 
sensors are inserted inside the beads rather than at the points of contact [see panel (b) of Fig. [T], we averaged the 
force between adjacent beads, F^ = {Fj + Fj+i)/2, as discussed in detail in Ref. [l5i |. 

Our numerical and experimental results greatly expand the previous work reported for chains of dimers composed 
of particles with similar elastic moduli and different masses 12]. For example, we include investigations of dimers 
composed of particles with moduli of different orders of magnitude. We also study dimers composed of different 
numbers of particles {Ni : 1 and 1 : chains), thereby varying the "thickness" of the heterogeneous layers. 



A. 1:1 Dimers 



We begin by discussing our results for 1:1 dimers. In Fig. [21 we show experimental and numerical results for 1 : 1 
dimers of steehPTFE [panels (a,b)] and steehrubber [panels (c,d)] particles. The steehPTFE chain consisted of 38 
beads, and the steehrubber chain had 19 beads; in each case, we used a steel bead as the striker. The dynamics indicate 
that the initial excited impulse develops into a solitary wave within the first 10 particles of the chain. Observe in the 
steehPTFE numerics the presence of small-period oscillations after the large solitary pulse. These arise as residual 
"radiation" emitted by the tail of the wave as it attempts to progress through the highly nonlinear lattice. This is 
reminiscent of radiative phenomena in different classes of nonlinear lattices, as discussed, for example, in Ref. [45| . 

Interestingly, even the chains composed of alternating steel and Neoprene rubber beads (inherently nonlinear elastic 
components) support the formation of solitary-like pulses. However, this configuration is highly dissipative and 
the solitary waves are consequently short-lived (although this feature is not incorporated in the theoretical model 
considered herein). As shown in Fig.[21^c,d), the initially very short pulse is quickly transformed into a much wider and 
slower pulse compared to the dimer composed of steel and PTFE particles. Additionally, the presence of dissipation 
in the steehrubber configuration tends to dampen the propagation of the pulses after the first 15 beads. Incorporating 
dissipative effects is a natural direction for the refinement of the model to be pursued in future studies. 

To illustrate the robustness of solitary-wave formation, we also studied dimer chains composed of other types 
of materials, including steehbronze (33 total particles), PTFE:glass (34 particles), and PTFE:nylon (33 particles) 
configurations. Two important properties of the pulses observed in these systems are their propagation speeds and 
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FIG. 2: (Color online) Force versus time response obtained from chains of dimers consisting of 1 stainless steel bead alternating 
with 1 PTFE (a,b) or 1 rubber (c,d) bead. Panels (a,c) show experimental results, and (b,d) show the corresponding numerical 
data. For both configurations, the initial velocity of the striker was 1.37 m/s on impact. The y-axis scale is 2 N per division 
in (a,b) and 20 N per division in (c,d). The numbered arrows point to the corresponding particles in the chain. In (a,b), the 
second curve (showing the results for particle 24) represents a PTFE bead and the other curves represent steel beads. In (c,d), 
particle 1 is steel and particle 6 is rubber. The force on the rubber particle is magnified by a factor of 20 for clarity. 



widths (measured by the full width at half maximum, or FWHM). We compute the pulse speed using time-of-flight 
measurements for both experiments and computations. To do this, we first measure the peak force of the pulse at 
each bead and determine the times at which these peaks occur. We then estimate the velocity of the pulse at each 
bead by examining how long it takes for the pulse (as measured at the peak time) to travel from an earlier bead to a 
later one. For the numerical calculations, we use the same bead separation distance as in corresponding experiments. 
For example, for the 1:1 steekPTFE dimer chain experiments, forces were measured at beads positions 13 and 24, so 
in the numerics we use time-of-flight measurements with an interval of 11 beads. 

We depict the results for steehbronze in Fig. [31[a,b) and for PTFE:glass in Fig. [3I^c,d). The steel-bronze dimer 
chain included 38 total particles; its striker was a steel bead, which was dropped from about 9.5 cm. We investigated 
this chain in order to examine a configuration with two different types of large-mass materials with relatively large 
elastic moduli. It is evident from both experiments and numerical simulations that the system supports the formation 
and propagation of highly nonlinear solitary waves. These waves have a pulse speed of 499.2 m/s experimentally 
and 544.0 m/s numerically. Their FWHM was measured at 2.03 beads and 2.06 beads in experiments and numerics, 
respectively. PTFE:glass dimer chains provide another configuration consisting of two different particles with similar 
masses, but in this case the masses are smaller than with steehbronze and the elastic moduli and (especially) Poisson 
ratios of the two materials are rather different. Although it is composed of lighter and softer particles, this dimer chain 
nevertheless supports the formation and propragation of solitary-like pulses. However, many of their properties are 
quite different: the signal speed is significantly reduced (it is 151.7 m/s in experiments and 221.4 m/s in the numerics) 
and the pulse FWHM is 1.99 experimentally and 2.10 beads numerically. In addition, secondary oscillations can be 
observed both in experiments and (to a lesser extent) in numerical simulations [see, in particular, the bottom (wall) 
signal in Fig.[2I^d)]. Finally, the PTFE:nylon chain (with a PTFE striker dropped from a height of 2.5 cm) provides a 
configuration with one small mass and a second mass that is extremely small. Unlike rubber, the nylon beads do not 
lead to severe dissipation. However, this configuration appears to have greater difficulty in forming and supporting 
solitary waves, at least in the chain lengths studied (33 beads). 

As we have suggested, dimer chains have decidedly different properties from uniform chains. We find numerically 
that the solitary pulse in uniform chains has an FWHM of about 2.06 cells (i.e., 2.06 beads), consistent with prior 
observations (see, e.g., Fig. 1.4 in ^1^]). Because of the similar masses in the steehbronze and PTFE:glass nylon chains. 
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FIG. 3; (Color online) Force versus time response obtained from chains of dimers consisting of (a,b) 1 steel bead alternating 
with 1 bronze bead and (c,d) 1 PTFE bead alternating with 1 glass bead. Panels (a,c) show experimental results, and (b,d) 
show the corresponding numerical data. The numbered arrows point to the corresponding particles in the chains. For the 
steehbronze configuration, a steel striker impacted the chain with an initial velocity of 1.21 m/s. The y-axis scale is 5N per 
divison in (a) and 20 N per division in (b). For the PTFE:glass configuration, a glass striker impacted the chain with an initial 
velocity of 1.17 m/s. The y-axis scale is 1 N per division in (c) and 3 N per division in (d). 



the pulse width in those 1:1 dimer chains is similar to that obtained for uniform chains. We observe this feature both 
experimentally and numerically. (As discussed later, the theoretical prediction for the mean FWHM in dimer chains 
depends only on the mass ratio of the two bead materials; the existence of the dimer is built into the long-wavelength 
asymptotics used to derive this result.) This suggests that there is a "critical" mass ratio 7711/7712 necessary for the 
system to become "sensitive" to the existence of the dimer. This poses an interesting question concerning what is the 
critical point (mass ratio) for the dynamics to resemble that of a chain with periodic "defects." That is, how large (or 
small) should 7712 be to obtain bona fide dimer dynamics (increased pulse width, etc.)? 

In the case of homogeneous chains, asymptotic analysis predicts a wave width of about 5 cells from tail to tail p^ . 
The propagation speed of the pulses in steel is quite large; for steel strikers with an impact velocity of 1.37 m/s, we 
found numerically that the pulse propagates at almost 700 m/s, and other authors have found large experimental 
propagation speeds with other initial conditions In contrast, as discussed in detail below, pulses in dimers 

have a shorter width in terms of number of cells. With a steel striker with a velocity of 1.37 m/s on impact, we 
find a propagation speed of 168.1 m/s numerically and 127.8 m/s experimentally. We discuss the width and pulse 
propagation speed of solitary waves in more detail below and include an asymptotic analysis of the width of 1:1 dimers 
as compared to the previously-studied limiting cases of monomer chains and 1:1 dimers with beads of mass mi and 
1712 ^ mi p^ . Preliminary experiments and numerical simulations on 1:1 dimers, considering only materials with 
elastic moduli of the same order of magnitude (and mass ratios 7711/7712 of 2, 4, 16, 24, and 64), were described in 
Ref. [H. 



B. TVi : 1 and 1 : iV2 Dimers 



We consider configurations with different periodicities by varying the number of steel and PTFE particles in steel: 
PTFE dimers. We obtain robust pulses for Ni : 1 dimers with TVi > 1 (with Ni as high as 7 in experiments and as 
high as 22 in numerics), though the transient dynamics and spatial widths of the developed solitary-like waves are 
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FIG. 4: (Color online) Force versus time response obtained from chains of dimers consisting of (a,b) 2 and (c,d) 5 stainless steel 
beads alternating with 1 PTFE bead. The total number of beads is 38 in each case. Panels (a,c) show experimental results, 
and (b,d) show the corresponding numerical data. For both experimental configurations, the striker had an impact velocity of 
1.37 m/s with the chain. The y-axis scale is 2 N per division for (a,c) and 1 N per division for (b,d). The numbered arrows 
point to the corresponding particles in the chain. For both configurations, the second curve (showing the results for particle 
24) represents a PTFE bead and the other curves represent steel beads. 



different. This is clearly illustrated in Fig. |4] through experimental and numerical results for 2 : 1 and 5 : 1 steehPTFE 
dimers. The difference in transient dynamics can be intuitively understood in terms of the increasing "frustration" of 
the system as Ni increases (due to the presence of the 1 bead of different type in each cell of an otherwise uniform 
chain). On the other hand, the spatial widths differ through a mechanism similar to the one that we will explain 
analytically (see the discussion below) for the 1 : 1 configurations. Each bead in a given cell responds differently but 
consistently throughout the chain. For example, in the 2:1 chain, the second steel bead in a cell possesses a cleaner 
and more solitary-like shape than the first or third bead. This is especially evident in the numerical simulations. 

For the 5:1 dimers, there is a critical "frustration" between the pulse size one would obtain in a uniform chain of 
beads (for which the predicted width is « 5 particles long J^]) versus that one would obtain in a 5:1 chain (about 30, 
corresponding to the presence of 5 unit cells with 6 particles each. The pulses are of smaller (spatial) width in terms of 
the number of cells than they are in the 1:1 dimers. More generally, we observed both numerically and experimentally 
that the width (in terms of the number of cells) of the solitary waves decreases gradually with increasing numbers A^i 
of steel beads in a single dimer (cell) from the 5-cell width expected for Hertzian interactions. It takes roughly six 
cells to achieve a stable pulse for the 2:1 configuration, whereas it takes 12 cells for the 5:1 configuration. (It also takes 
roughly 6 cells for the pulse in the 3:1 chain to stabilize and roughly 11 cells for that in the 4:1 chain to stabilize.) 
The increased propagation distance needed to stabilize the pulse is a result of the frustration discussed above. This 
frustration reflects the increased difficulties in considering an TVi : 1 cell as a "quasiparticle" as A^i increases. 

We observe in our experiments that the width of the solitary waves (in terms of the total number of cells) decreases, 
although it increases in terms of number of "participating" beads. For the numerics, the FWHM seems to be steady at 
first, but it is smaller for A^i = 4 as the "frustration" arises. We observe that the pulse speed increases with increasing 
A^i in both experiments and numerics. This is physically justified because the increase in the number of steel particles 
makes the system generally "stiffer." For each configuration, the scaling between the pulse's maximum force and 

its speed of propagation is Vs ~ ^l/^i whose validity we verify analytically for 1:1 dimers using long- wavelength 
asymptotics (see below). We compare the shapes of the experimental and numerical pulses at equilibrium for example 
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FIG. 5: (Color online) Rising shape comparison of two different solitary waves recorded in the middle of the chain for different 
dimer configurations. Each panel displays the evolution of the force (in N) with time (in microseconds). Curve 1 corresponds 
to the experimental points, and curve 2 represents the numerical values. The depicted examples are (a) steel-bronze dimer 
curves detected in particle 12 and (b) PTFE-giass dimer curves detected in particle 12. Note that in order to compare the 
shapes, the experimental curves have been multiplied by an arbitrary factor to obtain the same force amplitude as that found 
in the numerical data (where no dissipation is present). 

1:1 chains in Fig. [51 

We performed analogous experiments and numerics for relatively short chains consisting of 38 particles assembled 
in a 1 : A^2 setup (i.e., 1 stainless steel and N2 > 1 PTFE particles). We show the results for 7V2 = 2 in Fig. [ni^a,b) 
and those for A^2 = 5 in Fig. [6ljc,d). For N2 — 2, we again observed solitary waves, though some additional smaller 
amplitude structures were also present. For N2 = 5, we did not obtain robust, individual solitary pulses in either 
experiments or numerical simulations of short chains but instead observed a complex "rattling" between the particles, 
which resulted in the formation of trains of pulses of different amplitudes. Some structure was more apparent in the 
longer 1 : N2 chains we simulated numerically, but we still did not obtain robust solitary waves (see the discussion 
below) . 

To confirm that robust solitary waves form even for dimers with larger numbers of steel particles in a cell (i.e., large 
A^i in iVi : 1 steehPTFE dimer chains) despite the "frustration" phenomena observed above, we conducted additional 
numerical investigations with a lar ge n umber of particles. We used horizontal chains of beads (setting g = 0) to avoid 
the "linearizing" effects of gravity [13, HI] due to the significant gravitational precompression arising from the large 
number of particles. 

In panel (a) of Fig. [71 we show a space-time plot for a 15:1 steehPTFE chain of 1000 particles. As the figure 
demonstrates, one gets well-formed localized pulses even for larger iVi when the chain is sufficiently long. However, 
the shape of the pulse never stabilizes fully, as there are fluctuations in the peak as well as more extensive fluctuations 
next to the peak than what one sees near the peaks for smaller A^i (which do stabilize). Panel (b) of Fig. [71 shows 
a space-time plot for a 1:10 steehPTFE chain of 500 particles. While a larger peak is clearly observable, one does 
not obtain a pulse that is significantly larger than its surrounding peaks the way one does for dimers with several 
consecutive large-mass particles. We thus conclude that dimers with several consecutive large-mass beads behave 
rather differently in certain respects than dimers with several consecutive small-mass particles. This asymmetry 
would be interesting to study systematically in the future, as it might provide insights into considering the lattice 
chain discretely versus as a continuum. For example, we expect different excitation modes to be excited in lighter 
versus heavier particles. 
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FIG. 6: (Color online) Force versus time response obtained from chains of dimers consisting of 1 steel bead alternating with 
(a,b) 2 and (c,d) 5 PTFE beads. In each case, the total chain length was 38 particles. Panels (a,c) show experimental results, 
and (b,d) show the corresponding numerical data. For both configurations, the striker (a steel bead) had an impact velocity of 
1.37 m/s with the chain. The y-axis scale is 1 N per division in (a,c) and 2 N per division in (b,d). The numbered arrows point 
to the corresponding particles in the chain. In (a,c), both particle 14 and particle 27 are made of PTFE. In (b,d), particle 14 
is PTFE and particle 26 is steel. 




FIG. 7: (Color online) Space-time diagrams for long chains of (a) 15 : 1 and (b) 1 : 10 steel:PTFE dimers. For both 
configurations, the striker had a velocity of 1.37 m/s on impact. The striker was made of steel for the 15:1 dimer and made of 
PTFE for the 1:10 dimer. 
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C. Theoretical Analysis 



In our theoretical considerations, we investigate the prototypical dimer chain composed of 1 : 1 dimer cells containing 
beads with different masses (denoted by mi and 7712). We derive for the first time the long- wavelength equation for 
an arbitrary mass ratio and a general interaction exponent k. (We included a brief synopsis of this derivation in 



Ref. [30|.) The relevant rescaled dynamical equations can be expressed as [12 1 



miilj = {wj - Uj)'' - {uj - Wj^i)^ , (4) 
m2Wj = [uj+i - Wj)'' - {wj ~ Ujf- , (5) 

where Uj (wj) denotes the displacement of the jth bead of mass mi (1712). We present our analysis for the nonlinear 
wave in the case of general power-law interactions to illustrate the generality of the approach. We subsequently use 
Hertzian contacts (i.e., the special case of = 3/2) in order to compare our theoretical analysis with our experimental 
results. The fc = 1 case corresponds to linear interactions between particles. 

The distance from Uj to Wj (from Uj to uj^i) is denoted by D {2D). This gives a small parameter, allowing to 
us develop a long-wavelength approximation (LWA) by Taylor-expanding Eqs. (I4])-(I5]). In particular, we express Wj+i 
{wj-i) as a function of Uj (wj). The resulting PDEs (one for each mass "species") need to be "homogenized." To 
accomplish this, we follow 46], postulating a consistency condition between the two fields: 



w 



X{u + biDux + b2D^Uxx + hD^Uxxx + h^D^Uix + •••)■ (6) 



We then self-consistently determine the coefficients A and hi by demanding that Eqs. ([4]) and ([5]) be identical at each 
order. Note that the subscripts in the LWA denote derivatives. For some technical issues regarding the consistency 
condition, we defer the interested reader to the discussion in [i^ [see, in particular, Eqs. (3.2)-(3.7) and the accom- 
panying discussion]. The parameter A can take the values 1 (for acoustic excitations) or — mi/m2 (for optical ones). 
Given that our initial excitation generically produces in-phase (acoustic) waveforms, we restrict our considerations to 
A = 1 in the following exposition. Observe that for optical (out-of-phase) excitations to be produced, a fundamentally 
different and more complex experimental setup is necessary (in comparison to the experiments discussed herein) to 
be able to drive the system at a frequency within the gap of its linear spectrum. 

A direct comparison of the resulting PDEs for the two mass types at orders d'^-D'''^^ result in the following 
consistency constraints: 

61 = 1 , (7) 
h^^^, (8) 



mi -t- m2 
2mi — m2 
3(mi -I- m2) 



(9) 



mi{ml - mim2 + ml) 

04 = ; 7^5 ■ (10) 

3(toi -I- m2)'^ 

We note in passing that these results accurately capture the uniform-chain limit of mi = m2. In this case, Wj becomes 
effectively Uj+i and hence the consistency condition of ([5]) becomes the Taylor expansion oiuj. The quantity D, which 
equals the sum of the radii of consecutive beads in the case of dimer chains, is the radius of this expansion. 

The resulting PDE bearing the leading-order discreteness corrections (i.e., incorporating terms that are below 
the leading-order term) is of the form 

Urr = ul'^Uxx + Gul^^ul^ + Hu^^^UxxUxxx + lu^'^Uix , (H) 



where r = tyj2kD^+^ /{nii -t- m2) is a rescaled time. The constants arising in Eq. (fTT|) are given by 

^^^2(2^3fc + PK^ (12) 



^^j^2 - l)(2m^ + mima ~ ml) ^^^^ 



6(mi -f m2)2 
mim2 -|- 
6(mi -I- m2)^ 



T 2(m? — mimo -|- mn) 
^ = ^-^-4? To ■ (14) 
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Given the nature of the excitation propagating through the chain of beads, we seek travehng wave solutions u = u(^), 
where ^ = a; — Vst is the standard traveling wave variable with (renormalized) speed Vg — This ansatz results in 
an ordinary differential equation (ODE) for — v [this change of variable reduces the equation to third order]. In 
all the integrations that follow, we consider solutions with homogeneous Dirichlet boundary conditions far from the 
wave. 

We now need to use some tools from the theory of ODEs to further reduce the order of the equation. In particular, 
we change variables using v — , where the power p is chosen so that terms proportional to zP~'^z| disappear for 
z = z{i): 



{H + 3/) ± ^{H + 3/)2 - 8J(G + H + I) 
^'■' = 2(G + g + 7) • 

Observe that there are two possible values of p that lead to this reduction, but this does not change the essence of 
the results in what follows. We obtain an ODE of the form 

V'^pz^-^z^ = pz'^P-iz^ + p/z'^P-iz^j^ + [3p(p - 1)/ + p^ij] z'=P-22jZ{{ . (16) 

In Eq. (fTB]) . we use an integrating factor z", with 

a = l-fcp + 3(p-l)+piJ//, (17) 



to convert it to a tractable second order ODE, 



z^^ = /iz'' — az , (18) 



where \i = V]; + a)], ry = 1 + p — fcp, and a — l/[I{kp + a)]. One can find exact periodic solutions of Eq. (fT8|l 
that take the form 

= V = z" = B cos^ (/3C) , (19) 

where B = {^l/[(3'^s{s - /? = ^/^(l - 77)72, and s = 2/(1 - 77). 

The coexistence of such an exact trigonometric solution with nonlinear dispersion in the LWA of our model sug gest s 
the potential existence of essentially finite- width solutions that consist of a single arch of the profile of Eq. ^9]) [13 | . 
Two key experimentally testable properties of such solutions are as follows: 

• The amplitude- velocity scaling B ~ Vs^'"^ Observe that this scaling does not change from the single-species 
case [12, [l3|- Accordingly, we conclude that this is a geometrical result that arises from the power of the 
nonlinearity (i.e., from the geometry of the contacts between consecutive particles). 

• The solution width that is equal to 7r//3 and which, consequently, depends directly on the mass ratio. 
Among these two properties, the width is the one that naturally showcases the relevance and novelty of the general 



dimer theory developed herein in comparison with the monomer theory that can be found in IJ] . The advantage of 
the present theory is that it allows us to obtain an explicit, closed-form expression for the width that is valid for all k. 
Although the formula is too lengthy to write explicitly in the general case, we present here the expression for /? in the 
experimentally relevant caseof fc = 3/2: /3 = (30-8w + 8a;2)-i(V3(8-5w-h5u;2- V4 - 4w 13a;2 - ISw^ + 9w4)(((l4- 
uf{lb~'ii^ + 'ii^'^))/{D^{M-A2uj + muj^-Mi^^ + 17u^ + {-^ + buj-^M^)^/A-Au + 12,u'^ - 18^3 9w4)))i/2)^ ^j^^j.^ 
uj = milm\. One of the appealing features of this result is that it naturally generalizes all previously- known limiting 
cases - namely, the monomer case (mi = 7712), for which /3 = \/lO/(5Z3) [resulting in pulses extending to Stt/a/IO w 5 
sites], and the case with one mass that is much larger than the other (toi ^ 7712) [resulting in pulses of a/IOtt « 10 
sites, consisting of 5 cells with 2 sites each] [l^ . 

Before presenting a comparison between our theoretical and experimental results, we wish to briefly highlight that 
there is an inherent asymmetry in the quantity /3 in Eq. (|19p . This quantity is scaled to the spatial wavelength of 
D = Ri + R2 and depends intrinsically on the ratio of the radii and masses of consecutive particles through the 
parameter combination P2/P1 x ^2/^1 (where pj is the density of particle j). Consequently, one should not expect 
the resulting expression for (3 to be symmetric about lu ^ 1, so that it matters whether Ri or i?2 is larger. We show 
this asymmetry in Fig. [8] for the case of Hertzian interactions. Observe in the inset, however, an approximate local 
symmetry about uj-^^^ ~ 0.5 and the consequent consistency in the derived wavelengths of lOD w 10i?i at uj — 
(where (3 = l/-\/lO) and 5D — 5(i?i -I- R2) = 10i?i about lu = 1 (where /3 = l/-\/TO again). This also brings up the 
question of when it is appropriate to employ continuum versus discrete treatments of the granular chains. 
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FIG. 8: The parameter /3 as a function of the dimensionless mass-ratio parameter u^^^ for a system with k = 3/2, pij p\ — 1, 
i?i = 1, and i?2 £ [0, 10]. Note the approximate local symmetry in the inset. 

D. Comparison between experiments, numerical simulations, and theory 

First, we discuss a relation for the scaling of the maximum dynamic force Fm of the solitary waves versus their 
propagation speed Vs- From the theory, we find that Fm ~ ~ yj2fc/(fc-i)^ which for fc = 3/2 yields Vg ^ Fm^ , 
the same relation as for monomer chains [l2|. We tested this result both experimentally using 1:1 steehPTFE dimer 
chains and numerically for both 1:1 steehPTFE and 1:1 steehrubber chains. In all cases, the chains have 38 total 
beads. In the experiments, we dropped a striker (a stainless steel bead) from heights between m and 1.2 m and 
measured the peak force in beads 15 and 34. We averaged their amplitudes [Fm = (F15 + i^34)/2] and obtained 
the corresponding wave speed (Vs) using time-of- flight measurements. We computed Fm and Vs in exactly the same 
manner in our numerical simulations and obtained good agreement, as indicated in Fig. [9l^a). We show error bars for 
the experimental results. The errors in the force amplitude stem from oscilloscope noise and sensor orientation (their 
calibration varies a bit because they are aligned with the horizontal cut); those in the propagation speed arise from 
oscilloscope noise and the difficulty in determining exact time-of-flight distance between the peaks. In conducting 
our numerical simulations, we performed a least-squares fit to the obtained force and velocity values. To compare 
the numerical results with the theory, we repeated the numerical simulations in the absence of gravity, obtaining 
Vs ~ Fm^^^^ , which is in excellent agreement with the theory (even though we used a small number of particles). We 
remark that the inclusion of gravity, which decreases the exponent 7 in ~ F^rn ^ larger effect on the force- 
velocity scaling in dimer configurations that include very light particles (such as steehrubber dimers). We verified 
numerically that the same force-velocity scaling also holds for A^i : 1 steehPTFE dimer chains with larger iVi and for 
1:1:1 trimer chains. 

We saw earlier that the FWHM of solitary waves in a dimer chain is larger in terms of number of particles (though 
smaller in terms of number of cells) than that in a monomer chain. The physical mechanism behind this observation 
can be explained using the "homogenized" theory (that is, the long-wavelength asymptotics) derived above. We 
demonstrate this by examining the evolution of the FWHM as the wave propagates down the chain of beads. As 
shown in Fig. HKb) , the numerical and experimental results are in very good agreement with each other and with the 
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FIG. 9: (Color online) Comparison of experiments, numerical simulations, and theory, (a) Scaling of the maximum dynamic 
force Fm versus pulse propagation speed Vs in a 1:1 stainless steehPTFE dimer chain. The numerical simulations are shown 
by the solid curve and the experimental results are shown by points (with error bars). The dashed curve shows the numerical 
results using the elastic modulus E = 0.6 GPa for PTFE, the nominal static value reported in most of the literature [isl. Isoj]. 
(b) Evolution of solitary wave width (the FWHM) as a function of bead number. The experimental values are shown by solid 
(green) squares and the numerical values are shown by open (red) circles. (In both cases, we include curves between the points 
as visual guides.) The theoretical value for the FWHM with mi ^ m2 is given by line (1), that for the 1:1 steehPTFE chain 
is given by line (2), and that for a homogeneous chain is given by line (3). 



long-wavelength approximation that captures the average of the FWHM oscillations. The maximum dynamic forces 
(Ffn) in the numerics and experiments alternate from one bead to another because consecutive beads are composed 
of different materials. From Fig. [Hb), it is clear that the average FWHM of the dimer differs both from that of a 
homogeneous chain (mi = 1712) and from that of the mi 3> m2 limiting case discussed in [l^]. The average value of 
this width is best captured (for the proper ratio of masses) by the homogenized theory derived above, clearly evincing 
the relevance of our theoretical approach. To provide a quantitative measure of the theory's accuracy, we compute 
the relative error of the theoretical prediction versus the computational-average FWHM and compare it to those for 
the two limiting cases; we find errors of 22.8% for mi = m2, 9.2% for mi ^ 1112, and 2.5% for the correct mass ratio. 
The numerical simulations achieve a nearly-steady FWHM value by about the 4th bead, and the experiments achieve 
this by about the 6th bead, so the theoretical value for the width of the solitary wave becomes valid very early in the 



14 



Configuration 


Properties of Bead 1 


Properties of Bead 2 


Properties of Bead 3 


Steel:Bronze:PTFE 


large mass, large E, small v 


large mass, medium E, medium i/ 


small mass, small E, large v 


Steel:PTFE:Rubber 


large mass, large E, small i/ 


small mass, small E, large f 


very small mass, very small E, large f 


Steel:Glass:Nylon 


large mass, large E, small i/ 


small mass, medium E, very small v 


very small mass, small E, medium f 



TABLE II; Summary of the trimer configurations that we investigated in terms of their relative masses, elastic moduli E, and 
Poisson ratios i/. 



wave propagation dynamics. 

We now comment briefly on the relative importance of gravitational effects for chains of beads with small and large 
numbers of particles. First, because of the chain lengths and the pulse amplitudes considered in our experimental 
configurations, the pulse propagation can be safely assumed to be within the highly nonlinear regime. We consequently 
neglected the effects of gravity in the theoretical analysis, in accordance with the discussion of Ref. [2^. However, 
for vertical chains of beads that contain a large number of particles (several hundred or more) or are excited by 
smaller-amplitude pulses, the presence of the nonuniform gravitational precompression should be taken into account. 
As discussed in Ref. [l^ , this leads to a 1/3 power-law scaling in the wave width as a function of particle number as 
well as a power-law drift in the pulse speed. 



V. CHAINS OF TRIMERS 



To examine more intricate heterogeneous chains, we also performed experiments and numerical simulations for 
various trimer configurations. In Fig. llOf a). we depict the basic experimental setup of a trimer chain in which each 
cell consists of A^i consecutive beads of one type, N2 consecutive beads of a second type, and N3 consecutive beads of 
a third type. We focused on iVi — N2 = N3 = 1 and considered various material components with different masses, 
elastic moduli, and Poisson ratios. The 1:1:1 configurations we investigated were composed of steel:bronze:PTFE (steel 
striker, 36 particles), steel:glass:nylon (steel striker, 33 particles), and steel:PTFE:rubber (steel striker, 27 particles). 
Our specific choices of configurations were motivated by the desire to examine various combinations of hard, soft, 
heavy, and light materials (see the summary in Table HI)) . 
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FIG. 10: (Color online) (a) Experimental setup for trimer chain consisting of a periodic array of A''i consecutive beads of one 
material, consecutive beads of a second material, and consecutive beads of a third material. (b,c) Force versus time 
response obtained from chains of trimers consisting of 1 steel particle, 1 bronze particle, and 1 PTFE particle. Panel (b) shows 
experimental results and (c) shows numerical ones. A steel striker was dropped with a velocity of 0.54 m/s on impact, and the 
t/-axis scale is 1 N per division in both plots. The numbered arrows point to the corresponding particles in the chain. The 12th 
particle is PTFE, the 20th particle is bronze, and the 28th particle is steel. 

As with the dimer chains, we calculated the pulse speed for each configuration using time-of-flight measurements 
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FIG. 11: (Color online) Rising shape comparison of two different solitary waves recorded in the middle of the chain for different 
trimer configurations. Each plot displays the evolution of the force (in N) with time (in microseconds). Curve 1 corresponds 
to the experimental points, and curve 2 corresponds to the numerical ones. The depicted examples are (a) steel-glass-nylon 
trimer curves detected in particle 13 and (b) steel-bronze-PTFE curves detected in particle 12. Note that in order to compare 
the shapes, the experimental curves have been multiplied by an arbitrary factor to obtain the same force amplitude as that 
found in the numerical data (where no dissipation is present). 

between the sensors. We obtained very good qualitative and reasonable quantitative agreement (very good, given the 
lack of dissipation in the models) between the results of experiments and numerics. In Fig. fTOT b.c). we show the force 
versus time response for a chain of 1:1:1 steel:bronze:PTFE trimers. We observed the formation of robust solitary 
waves with a pulse propagation speed between the first and second sensor of about 115 m/s experimentally and 160 
m/s numerically. The FWHM is about 4.2 particles (1.4 cells) experimentally and 4.3 particles (1.4 cells) numerically. 
This illustrates the very good agreement between experiments and numerical simulations, as also indicated by the 
direct comparison of the pulses in Fig. 1111 We also obtained very good agreement between experiments and numerics in 
the other configurations. In Fig.[T21 we show the force versus time response for chains of steel: glass: nylon trimers (left 
panels) and steel:PTFE:rubber trimers (right panels). The steel: glass: nylon trimer chain had a steel striker dropped 
from 4 cm; its pulse speed was about 220 m/s experimentally and 286 m/s numerically. The FWHM was about 1.6 
cells both experimentally and numerically. The steel:PTFE:rubber trimer chain had a steel striker dropped from 15 
cm; its pulse velocity was about 33.9 m/s experimentally and 56.0 m/s numerically. (Recall that the dissipative effects 
are larger for rubber than for other materials.) The FWHM was about 1.4 cells both experimentally and numerically. 

For each configuration, we robustly demonstrated the formation and propagation of solitary waves. These waves can 
be clearly discerned both experimentally and numerically for the steel:bronze:PTFE and Steel:PTFE:nylon chains but 
only numerically for the trimer chain containing rubber, in which dissipation is much more important. The presence 
of small secondary pulses in the configurations can also typically be observed. Additionally, one can see in both the 
experiments and the numerical simulations that the pulses are not perfectly symmetric near the top of the chains. 
This skewness disappears as a pulse moves down a chain. 

VI. CONCLUSIONS 

We examined the propagation of solitary waves in heterogeneous, periodic granular media consisting of chains of 
beads using experiments, numerical simulations, and theory. Specifically, we investigated ID granular lattices of 
dimers and trimers composed of varying numbers of particles (yielding different periodicities) that were made of 
different materials. We found both experimentally and numerically that such heterogeneous systems robustly support 
the formation and propagation of highly localized nonlinear solitary waves, with widths and pulse propagation speeds 
that depend on the periodicity of the chain (that is, the length of its unit cell). 
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FIG. 12: (Color online) Force versus time response for chains of 1:1:1 trimers consisting of (a,b) steel:glass:nylon and (c,d) 
steel:PTFE:rubber. Panels (a,c) show experimental results and panels (b,d) show numerical results. In (a,b), a steel striker 
was dropped with a velocity of 0.89 m/s on impact; the y-axis scale is 1 N per division in (a) and 2 N per division in (b). In 
(c,d), a steel striker was dropped with a velocity of 1.72 m/s on impact; the y-axis scale is 0.5 N per division in both (c) and 
(d). The numbered arrows point to the corresponding particles in the chain. In panels (a,b), both depicted beads are made of 
steel; in panels (c,d), the 7th bead is made of steel and the 17th is made of PTFE. 

For dimer chains consisting of cells composed of Ni particles of massive materials such as steel doped by 1 light 
particle such as PTFE, we find that the width (expressed as the number of unit cells) decreases, whereas the propa- 
gation speed increases with A^i. We also observe a "frustration" phenomenon for Ni > 4, although robust localized 
pulses nevertheless form even in this case. We showed numerically, experimentally, and analytically (using a long- 
wavelength approximation) that the force-velocity scaling for chains of dimers follows the same relation as that for 
homogeneous chains. We also found very good agreement for the solitary-wave widths (which depend on the mass 
ratio of the dimer materials) using all three approaches, generalizing previous results from all known limiting cases. 
Finally, we demonstrated both experimentally and numerically the formation and propagation of robust, strongly 
localized solitary waves in trimer lattices consisting of chains of cells composed of 1 particle of each of three varieties. 

The qualitative and quantitative understanding of the dynamics of dimer and trimer chains obtained in the present 
work, and their characterization as possessing highly nonlinear solitary waves with similarities (e.g., force- velocity 
scaling) as well as differences (e.g., in the width) from the standard monomer chain, paves the way for studies in 
increasingly heterogeneous media in both one and in higher dimensions. Such studies are currently in progress and 
will be reported in future publications. 
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